Systematics of g factors of 2+ states in even-even nuclei from Gd to Pt: 
A microscopic description by the projected shell model 



Bao-An Bian 1,2 , Yao-Min Di 1 , Gui-Lu Long 3,1 , Yang Sun 4 ' 1 , Jing-ye Zhang 5 , Javid A. Sheikh 6 

1 Department of Physics, Xuzhou Normal University, Xuzhou, Jiangsu 221009, P.R. China 
2 Institute of Low Energy Nuclear Physics, Beijing Normal University, Beijing 100875, P. R. China 
3 Department of Physics, Tsinghua University, Beijing 100084, P R- China 
4 Department of Physics and Joint Institute for Nuclear Astrophysics, 
University of Notre Dame, Notre Dame, Indiana 46556, USA 
5 Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 
6 Department of Physics, University of Kashmir, Hazrathbal, Srinagar, Kashmir 190 006, India 

(Dated: February 9, 2008) 

The systematics of g factor of first excited 2 + state vs neutron number TV is studied by the 
projected shell model. The study covers the even-even nuclei of all isotopic chains from Gd to Pt. g 
factors are calculated by using the many-body wavefunctions that reproduces well the energy levels 
■») . and B(E2)s of the ground-state bands. For Gd to W isotopes the characteristic feature of the g 

f~>) ■ factor data along an isotopic chain is described by the present model. Deficiency of the model in 

' the g factor description for the heavier Os and Pt isotopes is discussed. 

Ch . PACS numbers: 21.60.Cs, 23.20.-g, 27.70.+q, 27.80. +w 
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INTRODUCTION 



Nuclear magnetic dipole moment can provide valuable information on the microscopic structure of a nuclear system. 
, It is a sensitive probe of nuclear wavefunctions and hence can serve as a strict testing ground for theoretical models. 
Because of the intrinsically opposite signs of the neutron and proton g Sl a study of the gyromagnetic factor (g factor) 
enables determination of the detailed structure for underlying states. For example, variation of g factors often is a clear 
indicator for a single-particle component that strongly influences the total wavefunction. With advances in modern 
experimental techniques and sensitive detectors, progress in g factor measurement has continuously been made. In 
a recent paper [lj], Berant et at, by summarizing their new and the accumulated data from Refs. 0, Q, raised an 
' interesting question on the systematic behavior of the first excited 2 + g factors (denoted as g(2±) hereafter) in the 
rare earth nuclei and the heavier mass region. The data for these even-even nuclei indicate characteristic features of 
the systematics (see Fig. 4 below): with increasing neutron number, g(2f) factors display a decreasing trend in the 
Gd, Dy, and Er isotopes; stay nearly constant within a range of the Yb and Hf chains; then change to an increasing 
trend in the W and Os isotopes; but show a flat behavior in the Pt chain. 

Clearly, the overall behavior of these g factors exhibits a large deviation from the rotor value, Z/A, which has 
only a very weak and smooth dependence on nucleon numbers [4|. On the other hand, the proton- neutron version 
of the Interacting Boson Model [a] gives an overly strong particle number dependence, and thus fails to reproduce 
the fiat behavior of the g factors near the midshell [l[. These facts may suggest that in realistic nuclear systems, g 
factors reflect a delicate interplay between collective and single-particle degrees of freedom, which is dictated by the 
detailed shell structure. Very recently, Zhang et al. have discussed the systematical behavior of these g(2^) factors in 
terms of a simple phenomenological model [6(. However, it is desired that the g(2f) systematics can be described by 
microscopic theories. Spherical shell model calculation is applicable only to those heavy nuclei near the shell closures 



(for a recent example, 
heavy, deformed nuclei 



been microscopic models employed in the g factor calculation for 
However, except in Ref. [ll(, those calculations focused mainly 



on one or a few chosen examples in an isotopic chain. A microscopic description of the experimental g(2^) data for 
the large mass region as presented in Ref. [lj remains as a challenge to microscopic theories. 

In the present article, we carry out a systematical study for g(2f) factors. As a theoretical tool, we employ the 
code developed in Ref. [12j, which is based on the Projected Shell Model (PSM) [l6| . In the PSM approach, one 
introduces axially deformed basis and performs exactly angular-momentum projection on the intrinsic states from 
deformed mean-field calculations. For even-even nuclei, angular-momentum projection on the lowest K = state 
generates a rotational band, which is the main component of the ground-state band (g-band) including the 2^ state of 
our interest. The lowest K — state is the quasiparticle vacuum obtained microscopically through the Nilsson + BCS 
calculations in a large single-particle space. We thus expect this model to be an appropriate microscopic theory for 
the present investigation. The calculation is performed for nuclei from the Gd (Z = 64) to the Pt (Z = 78) isotopic 
chains, with neutron numbers ranging from N — 88 to 120. The nuclei studied here are known to have very different 
collective properties; for example, they can be well-deformed, less-deformed, or soft nuclei. The g factor of the first 
2 + state is dependent on details of the total wavefunction. The deviations from the collective Z/A trend are mainly 
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understood as a consequence of single-particle make-up of the wavefunction and the interplay between the collective 
degree of freedom and single particles. The present model employs deformed Nilsson single-particle states at fixed 
deformation. As we shall discuss, this scheme works well for well-deformed nuclei but is insufficient for description of 
soft nuclei in the heavier Os and Pt region. 

II. OUTLINE OF THE THEORY 

A calculation for medium to heavy deformed nuclei in terms of the conventional (spherical) shell model is not 
feasible despite recent computational advances. The successful description of heavy deformed nuclei can be traced 
back to the introduction of the Nilsson potential [13] • In the Nilsson model, nuclear states are described by considering 
nucleons moving in a deformed potential. Deformed states are defined in the body-fixed frame of reference in which the 
rotational symmetry is broken. In order to calculate the observables, it is necessary to restore the broken rotational 
symmetry in the wavefunction. This can be done by using the standard angular momentum projection method. The 
projected states are then used to diagonalize a two-body shell model Hamiltonian. Thus, our approach follows closely 
with the basic philosophy of the conventional shell model. The main difference to the conventional shell model is that 
in the PSM, one starts with a deformed basis rather than a spherical one. For the details of the projection method, 
we refer to the PSM review article [l6| and references cited therein. 

The PSM constructs the shell-model space by using the axially symmetric Nilsson states with a quadrupole defor- 
mation €2- Pairing correlations are incorporated into the Nilsson states by the BCS calculations. The consequence 
of the Nilsson-BCS calculations defines a quasiparticle vacuum |$(e2)) = \&) and the associated set of quasiparticle 
states in the intrinsic frame. For the low-lying nuclear states near the ground state the PSM wavefunction can be 
expressed as 

K I ) = f I K=oPL K =om, (i) 

which is the angular-momentum-projected quasiparticle vacuum state with / being the normalization factor. In Eq. 
(fTj), P[f K is the angular-momentum-projection operator [l8j 

Hik = J dsi D{ IK (n) R(Q). (2) 

As in the early PSM calculations, we use the pairing plus quadrupole-quadrupole (QQ) Hamiltonian with inclusion 
of the quadrupole-pairing term 

In Eq. ©, Hq is the spherical single-particle Hamiltonian, which contains a proper spin-orbit force [Ll\ . As dis- 
cussed in Ref. [l6|j the QQ- force strength x is adjusted such that the quadrupole deformation £2 is obtained as a 
result of the self-consistent mean-field HFB calculation. The monopole pairing strength Gm is of the standard form 
Gm = [20.12 =p 13.13(7V — Z)/A] /A, with "— " for neutrons and "+" for protons, which approximately reproduces the 
observed odd-even mass differences in this mass region [l6j]. The quadrupole pairing strength Gq is assumed to be 
proportional to Gm, with the proportionality constant being fixed to be 0.16 for all nuclei considered in this paper. 
The same constant has been used in the previous calculations for rare earth nuclei [l6[ • For the valence single-particle 
space, we have included three major shells, N — 4, 5, 6 (3,4,5), for neutrons (protons). 

In short, the procedure of the present calculation is that based on a deformed Nilsson potential with pairing 
included in the BCS treatment, one performs explicit angular-momentum projection with a two-body interaction 
which conforms (through self-consistent conditions) with the mean-field Nilsson potential. The Hamiltonian with 
separable forces serves as an effective interaction, the strengths of which have been fitted to experimental data. The 
deformed single-particle states with deformation parameters £2 are used solely as a starting basis. It is sufficient 
for a calculation to have these deformation parameters close to the "true" nuclear deformation. Of course, a large 
departure from a true deformation would result in a significant enhancement in dimension of the configuration space, 
with the extreme case being the conventional shell model based on a spherical basis (£2 = 0). 

III. ENERGY LEVELS AND BE(2) VALUES 

We study 61 nuclei from the Gd, Dy, Er, Yb, Hf, W, Os, and Pt isotopic chains. This large group includes nuclei 
with very different collective behavior. It is well-known that with neutron number around 90, nuclei are traditionally 
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TABLE I: The quadrupole deformation parameters with which the deformed bases are constructed. 



z 


= 64 


152 Gd 
0.212 


154 Gd 
0.278 


156 Gd 
0.295 


158 Gd 
0.305 


160 Gd 
0.320 












z 


= 66 

£2 


154 By 
0.200 


156 Dy 
0.240 


158 Dy 
0.260 


160 Dy 
0.270 


162 Dy 
0.280 


164 Dy 
0.290 










z 


= 68 

£2 


156 Er 
0.195 


i58 Er 
0.230 


160 Er 

0.257 


i62 Er 
0.265 


i64 Er 
0.258 


166 Er 
0.262 


168 Er 
0.273 


170 Er 
0.276 






z 


= 70 


164 Yb 
0.245 


166 

0.250 


168 Yb 
0.260 


i70 Yb 
0.265 


172 Yb 

0.269 


174y b 

0.266 


176 Yb 

0.263 








z 


= 72 

£2 


166 Hf 

0.219 


168 Hf 

0.235 


170 Hf 

0.240 


172 Hf 

0.250 


174 Hf 

0.253 


176 Hf 
0.246 


178 ffi 

0.241 


180 Hf 

0.231 






z 


= 74 

£2 


168 W 

0.193 


170 W 
0.201 


172 W 
0.217 


174 W 
0.220 


176 W 
0.225 


178 W 
0.195 


iso w 
0.190 


182 W 
0.195 


0.195 


186 W 

0.190 


z 


= 76 

£2 


178 Os 
0.188 


180 Os 
0.172 


182 0s 
0.170 


184 Os 
0.173 


186 Os 
0.158 


188 Os 
0.154 


190 Os 
0.150 


192 0s 
0.145 






z 


= 78 

£2 


182p t 

0.197 


184p t 

0.187 


186 p t 

0.175 


188 p t 

0.135 


190p t 

0.128 


192p t 

0.116 


194p t 

0.113 


196p t 

0.120 


198p t 

0.170 





4 



CD 



152 Gd D 
□ • 
• 

i 

1 * , , 


154 Dy 

i ' 

■ 

? 


166 Er 

m 

. . ' ' 


172 Yb 

■ 

».f, 


17S Hf 

■ 

. « ' 

1 — W — i — 1— — " — i 1 1 1 1 


180 w 

■ 

■ 

» f ' , ! 


186 0s 

■ 

■ 

■ 

. ■ 


n 
• 

i 

■ 

■ 



. » 

2 4 6 8 10 2 4 6 8 10 12 
Spin 



FIG. 1: Comparison of calculated energies with experimental data for g-bands. Open squares represent the calculated results 
and solid circles the data. 



known as 7-soft nuclei. On the other hand, the heavier isotopes in the Os and Pt chains contain also significant 
7-softness, and these are typical examples of 0(6) nuclei according to the Interacting Boson Model [l9j]. Between 
these two regions, nuclei are strongly deformed, for most of which the deformation is axial and the low-lying spectrum 
typically exhibits a rotor behavior. In Table I, we list the deformation parameters with which the deformed bases are 
constructed. The listed quadrupole deformations ti agree with the systematic trend of those experimentally adopted 
ones 20] although the absolute values of ours are smaller. Note that it is not necessary for our input deformations 
to be exactly the same as those extracted from experiment as long as the so-constructed bases can correctly describe 
experimental B(E2)s (see Figs. 2 and 3 below). In the calculation, when the calculation condition is fixed we do not 
have a freedom to readjust the parameters to reproduce the g-factors. Namely, under a fixed calculation condition, g 
factors are predicted and the underlying physics is discussed. 

Fig. 1 shows the calculated energy levels for g-bands together with experimental data. We include in the figure 
only one nucleus selected from each of the isotopic chains because all the calculations have achieved the same level 
of agreement. These examples are chosen to represent nuclei with distinct collective behavior. For instance, with 
neutron number 88, 152 Gd and 154 Dy are typical 7-soft nuclei lying in the transitional region. 166 Er, 172 Yb, and 
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FIG. 2: Comparison of calculated g-band B(E2) values with available data. Open squares represent the calculated results and 
solid circles (with error bars) the data. 



Hf are representative examples of strongly deformed nuclei lying in the midshell, which have nearly constant 
moment of inertia. Finally, 186 0s and 192 Pt are again 7-soft in the transitional region. For all these nuclei with very 
different rotational behavior, it can be seen from Fig. 1 that the energy levels have been well reproduced by the 
calculation. Note that for the strongly deformed nuclei lying in the midshell, the energy intervals are considerably 
smaller, corresponding to larger moments of inertia in these nuclei. The deviation seen in the higher-spin states in 
152 Gd, 154 Dy, and 192 Pt is caused by the g-band interactions with other configurations that have not been considered 
in Eq. p. 

We further calculate the g-band B(E2) values for all nuclei studied in this paper using the same deformation 
parameters listed in Table I. The B(E2) values that are related to an electric transition probability from an initial 
state J to a final state 1 — 2 are given by 

B(E2, 1 -> I - 2) = |( tt'- 2 ||Q 2 ||*' )| 2 , (4) 

where wavefunctions l^ 1 ) are those in Eq. (p}. The effective charges used in the calculation are the standard ones 
e w = 1.5e and e v = 0.5e. The effective charges are fixed for all nuclei studied in this paper without any individual 
adjustment. Any variations in the calculated B(E2) values, among one rotational band or between those in different 
nuclei, are subject to the structure change in wavefunctions. The calculation is compared with available data in Fig. 
2. Again, we include in this figure only one nucleus selected from each of the isotopic chains, the same set of nuclei as 
shown in Fig. 1. It can be seen that the B{E2) values are also nicely reproduced. Not only are the absolute values 
in each nucleus correctly given, but also the variations as a function of angular-momentum are described. We note in 
particular that the B(E2)s in 192 Pt are well reproduced; as we shall see later, for the Pt isotopic chain we encounter 
difficulties in the g-factor calculation. 

In Fig. 3, we plot the calculated B(E2, I = 2 — ► / = 0) values for all the 61 nuclei considered in this work, and 
compare them with available data. The numbers used for the figure are listed in Table II. We stress that the variations 
in B{E2) along each isotopic chain, i.e. a rapid increase up to N s=s 94, the flat behavior for 96 < N < 108, and a 
decreasing trend after the midshell, are correctly described. A few local exceptions with rather large B(E2) values in 
the data (such as in 172 W, 182 W, and 182 0s) can not be understood by the present calculation. The global trend of 
the B(E2) values with exclusion of the Os and Pt chains has been discussed by a simple one-parameter model [6j. 

The agreement of the calculated B(E2, I = 2 — > I = Q) values with data (Figs. 2 and 3) as well as the reproduction 
of the g-band energies (Fig. 1) indicate that under the present calculation conditions, we are able to describe the basic 
structure quantities for these nuclei. As far as the low-lying energy levels and B(E2)s are concerned, the method 
works well. The systematical agreement between the calculated and experimental B(E2)s (Fig. 3) justifies the use of 
the deformation parameters in Table I. 
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FIG. 3: Comparison of calculated B(E2, 1 = 2-^1 — 0) with available data. Open squares represent the calculated results 
and solid circles (with error bars) the data. 



IV. G-FACTOR RESULTS AND DISCUSSION 



We now turn our discussion to g factors. In the PSM, g factors can be directly computed as 



9(1) 



1 



li N I pt, N I 

with /J, T (I) being the magnetic moment of a state 



[/!*(/) + (I)] 



(5) 



Mr (7) = <*f|AII*f) 
I 



i 

7W+^ 



(^11^11*0 

gl^'WfW* 1 ) + (g T s - glWWW* 1 ) 



(6) 



where r = 7r and v for protons and neutrons, respectively. 

We use the same wavefunctions that are used to evaluate B(E2) values. In the angular-momentum-projection 
theory, the reduced matrix element for an operator m (with m being either j or s in Eq. can be explicitly 
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TABLE II: Comparison of calculated B(E2, J = 2 — ► I = 0) (in e 2 b 2 ) with available data. 



Z = 64 
Exp. 
Th. 


152 Gd 
0.35(3) 
0.57 


154 Gd 
0.76(2) 
0.92 


156 Gd 
0.93(3) 
1.03 


158 Gd 

1.01(3) 
1.09 


160 Qd 

1.03(3) 
1.18 






Z = 66 
Exp. 
Th. 


154 Dy 
0.48(4) 
0.55 


156 Dy 
0.75(1) 
0.77 


158 Dy 
0.93(4) 
0.90 


160 Dy 
1.00(3) 
0.97 


162 Dy 164 Dy 

1.04(3) 1.11(2) 
1.04 1.10 






Z = 68 
Exp. 
Th. 


156 Er 
0.33(1) 
0.52 


158 Er 

0.60(3) 
0.72 


160 Er 

0.86(4) 
0.89 


162 Er 

1.00(5) 
0.96 


164 Er 166 Er 168 Er 170 Er 

1.16(8) 1.16(5) 1.14(3) 1.16(2) 
0.97 1.02 1.10 1.14 






Z = 70 
Exp. 
Th. 


164 Yb 

0.93(3) 
0.83 


166 Yb 

1.03(5) 
0.90 


168 Yb 

1.15(4) 
0.96 


170 Yb 

1.12(3) 
1.03 


172 Yb 174 Yb 176 Yb 

1.20(1) 1.16(4) 1.07(4) 
1.08 1.10 1.07 






Z = 72 
Exp. 
Th. 


166 Hf 

0.69(4) 
0.66 


168 Hf 

0.84(4) 
0.77 


170 Hf 

1.01(30) 
0.83 


172 Hf 

0.88(6) 
0.89 


174 Hf 176 Hf 178 Hf 180 Hf 

0.88(5) 1.07(5) 0.95(3) 0.94(3) 
0.94 0.92 0.90 0.85 






Z = 74 
Exp. 
Th. 


168 W 

0.64(3) 
0.51 


170 W 
0.69(2) 
0.57 


172 W 
0.97(9) 
0.66 


174 W 
0.78(5) 
0.69 


!«W 178 w 180 w 182 w 184 w 186 w 

0.83(4) 1.06(2) 0.74(2) 0.70(1) 
0.73 0.65 0.63 0.64 0.62 0.59 


Z = 76 
Exp. 
Th. 


178 0s 
0.54 


180 Os 
0.72(18) 
0.51 


182 0s 
0.77(2) 
0.51 


1S4 Os 186 Os 188 Os 190 Os 192 Os 
0.61(15) 0.58(2) 0.51(1) 0.47(1) 0.42(1) 
0.49 0.45 0.43 0.40 0.37 






Z = 78 
Exp. 
Th. 


182p t 

0.51 


184p t 

0.70(3) 
0.45 


186p t 

0.59(3) 
0.39 


188p t 

0.53(9) 
0.26 


190p t 192p t 194p t 196p t 

0.36(2) 0.38(1) 0.33(1) 0.27(1) 
0.25 0.22 0.20 0.20 


198p t 

0.29 





expressed as 

<¥'||ih-||*') = E E {-Y~ Mi ( J M mm) (*\PK,M f miMP}c tMi \*) 

Ki,K f Mi ,Mf ,M V ' ' 

= (21+1) E I h fLfL. E 

Ki,K s M',M" 

x J dnD M »KM(®\miM>R(fy\®)- 
In our calculation, the following standard values for gi and g s appearing in Eq. ^ are taken: 

9l =1, 9s = 5-586 x 0.75, 
9i =0, g v 8 = -3.826 x 0.75. 

g™ and g v s are damped by a usual 0.75 factor from the free-nucleon values to account for the core-polarization and 
meson-exchange current corrections [2l|. The same values are used for all g factor calculations in the present paper, 
as in the previous projected shell model calculations, without any adjustment for individual nuclei. 

We present the systematics of g factor of first 2 + state for all the isotopic chains from Gd to Pt. In Fig. 4, a 
comparison of calculated results with available data [3, [3, Q is given. The numbers used for the figure are listed in 
Table III. Interesting systematical features are clearly observed. The experimental g{2f) values of Gd, Dy, and Er 
isotopes show a downsloping trend with increasing neutron number N; those of heavier W and Os isotopes exhibit an 
upsloping behavior; and, according to the current set of data, the g(2f) factors of Yb, Hf and Pt isotopes are almost 
constant within each isotopic chain. From Fig. 4, it can be seen that the observed trends for the isotopic chains 
from Gd to W are qualitatively reproduced by the present calculation. Especially for most Dy, Er, and W nuclei, a 
quantitative agreement with data is achieved. The downsloping trend of Gd, Dy, and Er isotopes is well described. 
The calculation further predicts that for the Yb, Hf, and W isotopic chains, a downsloping trend is still visible for 
N < 100, but becomes nearly constant at the neutron midshcll. The upsloping trend of W isotopes with N > 106 is 
well reproduced by the calculation. Nevertheless, there are cases where we find disagreement. We shall comment on 
those cases later. 
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FIG. 4: Comparison of calculated 2 + state g factor with available data. Open squares represent the calculated results and 
solid circles (with error bars) the data. 



Fig. 5 shows the calculated g(2f) factors for Er and W isotopes. To understand the iV-dependent variations, 
theoretical results are decomposed further into the individual contributions from the proton and neutron operator 
(see Eq. ((SJ)), the sum of which is the total g factor that was compared with data in Fig. 4. Remarkably, we find 
that the g(2^ ) variations originate mostly from the proton contribution while the neutron contribution is very small 
and stays nearly constant within each of the isotopic chain. Therefore, the change of the proton contribution alone 
describes the observed variation trends of g(2^) factors, namely, a rapid decrease followed by a constancy in the Er 
nuclei and probably a U-shape in the W nuclei (see Fig. 5). It is rather interesting that the proton contribution 
varies as the neutron number changes along an isotopic chain. This can only be possible with the presence of strong 
neutron-proton interaction. In this regard, we notice that Zhang et al. [6| suggested that the observed constancy of 
the g(2f) factors in the well-deformed region is attributed to the reduction of proton-neutron interaction strengths 
near the midshell. 

While the calculation predicts a flat behavior of g(2^) for the lighter Os nuclei, a clear departure from data occurs 
for heavier Os isotopes and for all the Pt isotopes considered in this paper. Although for 198 Pt, the theoretical value 
becomes closer to the data, we must conclude that the present calculation fails to describe the observed g(2± ) trend 
in the Os and Pt chains. We have tried various calculations by constructing our deformed model space with different 
£2 deformations. In the example of 192 Pt, it is found that with artificially increasing basis deformation, the g{2^) 
factor values starts going up, and at £2 = 0.24 the calculated g(2f) agrees with data. However, the experimental 
energy levels and B(E2) values in this nucleus cannot be simultaneously described. This may indicate that, although 
g-band energies and B(E2)s in the Os and Pt chains are reproduced by the model, the obtained wavefunctions with 
respect to the single-particle content can be wrong. We note that energy levels and B(E2)s near the ground state 
reflect mainly the collective properties of even-even nuclei and are not sensitive to single particles. 
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TABLE III: Comparison of calculated 2 + state g factor with available data. 



Z = 64 
Exp. 
Th. 


152 Gd 
0.48(4) 
0.43 


154 Gd 
0.48(3) 
0.31 


156 Gd 
0.41(7) 
0.29 


158 Gd 
0.39(3) 
0.27 


160 Gd 
0.36(2) 
0.26 












Z = 66 
Exp. 
Th. 


154 Dy 
0.36(2) 
0.47 


156 Dy 
0.39(4) 
0.39 


158 Dy 
0.36(3) 
0.34 


160 Dy 
0.35(2) 
0.31 


162 Dy 
0.35(2) 
0.28 


164 Dy 
0.34(1) 
0.26 










Z = 68 
Exp. 
Th. 


i56 Er 

0.40(6) 
0.49 


158 Er 

0.36(6) 
0.41 


160 Er 

0.33(6) 
0.35 


162 Er 
0.32 


i64 Er 

0.349(8) 
0.29 


166 Er 168 Er 170 Er 

0.325(5) 0.31(3) 0.317(7) 
0.25 0.27 0.26 






Z = 70 
Exp. 
Th. 


164 Yb 

0.32(3) 
0.32 


i66 Yb 
0.29 


i68 Yb 
0.25 


170 Yb 172 Yb 

0.337(4) 0.335(8) 
0.26 0.25 


174 Yb 176 Yb 

0.338(4) 0.34(2) 
0.23 0.25 








Z = 72 
Exp. 
Th. 


166 Hf 

0.32 


168 jj f 

0.28 


170 Hf 

0.24 


m Hf 
0.25 


174 Hf 

0.24 


176 ffi 

0.32(3) 
0.21 


178 Hf 

0.24(2) 
0.22 


IBOjjf 

0.31(2) 
0.24 






Z = 74 
Exp. 
Th. 


0.25(5) 
0.32 


«°w 

0.27 


172 W 
0.23 


174 W 
0.22 


176 W 
0.21 


178 w 

0.25(5) 
0.21 


180 W 
0.26(2) 
0.20 


182 w 

0.26(1) 
0.24 


184 W 186 W 

0.289(7) 0.31(2) 
0.24 0.26 


Z = 76 
Exp. 
Th. 


178 Os 
0.26 


180 Os 
0.26 


182 Os 
0.25 


184 Os 
0.25 


186 Os 
0.26(2) 
0.25 


188 Os 
0.29(1) 
0.26 


190 Os 
0.35(1) 
0.26 


192 Os 
0.40(1) 
0.26 






Z = 78 
Exp. 
Th. 


182p t 

0.11 


184p t 

0.28(3) 
0.11 


186 p t 

0.27(3) 
0.11 


188p t 

0.29(4) 
0.11 


190p t 

0.29(2) 
0.11 


192p t 

0.29(2) 
0.12 


194p t 

0.30(2) 
0.13 


196p t 

0.30(3) 
0.16 


198p t 

0.32(1) 
0.23 





This calls for a further improvement of the projected shell model type approaches to generally describe transitional 
nuclei. More correlations in the wavefunction need to be included, which goes beyond what an axially deformed 
quasiparticle vacuum state can contain. These correlations can be introduced by the addition of the D-pair operators 
to the vacuum state [23 |. which takes both quasiparticle and collective degrees of freedom explicitly into account in 
a shell model basis. The generator coordinate method, which consists of a construction of a linear superposition of 
different product wave-functions, can also be adopted. 

V. SUMMARY 

Inspired by the recent experimental work of Berant et at, [l| , we have made an attempt to study systematically g(2^) 
for all the isotopes from Gd to Pt, using the projected shell model approach. With a single set of interaction strengths, 
we have carried out calculations for each nucleus in a projected basis constructed with appropriate deformations. We 
have been able to reproduce the energy levels and B(E2)s for low-lying states in the g-band for all the 61 nuclei 
considered in the paper. With the same set of calculation conditions, we have calculated g factors of the first 2 + 
state. It has been found that for the isotopes from the Gd to W chain, the characteristics of experimental data along 
each isotopic chain are described by the theoretical calculations, such as the downsloping trend in the Gd, Dy, and 
Er isotopes, the upsloping trend in the W isotopes, and the flat behavior of the Yb and Hf isotopes. For the heavier 
Os and Pt nuclei, the results have indicated that although the energy levels and B(E2)s can be described equally 
well as in the lighter nuclei, the calculated g factors are wrong. Study of the separate contributions of proton and 
neutron to the g{2^) factors suggests that the variations of the g factors as the neutron number changes originate 
mainly from the proton contribution. The overly weak proton contribution to the g factors for Os and Pt isotopes 
indicates deficiency in the wavefunctions. To describe the heavier isotopes in the Os and Pt chains, improvement in 
the theory is required. 
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FIG. 5: Calculated first excited 2 + state g factors for Er and W isotopes, with decomposition of the total g factor (open 
squares) into proton (open circles) and neutron (open triangles) contributions. 
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